Cond-mat/0211466 



Theoretical analysis of the focusing of acoustic waves by two-dimensional sonic crystals 



Bikash C. Gupta and Zhen Ye[] 
Wave Phenomena Laboratory, Department of Physics, 
National Central University, Chungli, Taiwan 32054, R- O. 
(Dated: February 1, 2008) 



a 



(N 
O 
O 
(N 

> 
O 

(N 



> 



(N 
O 



C 

O 

o 



13 



Motivated by a recent experiment on acoustic lenses, we perform numerical calculations based 
on a multiple scattering technique to investigate the focusing of acoustic waves with sonic crystals 
formed by rigid cylinders in air. The focusing effects for crystals of various shapes are examined. 
The dependance of the focusing length on the filling factor is also studied. It is observed that 
both the shape and filling factor play a crucial role in controlling the focusing. Furthermore, the 
robustness of the focusing against disorders is studied. The results show that the sensitivity of the 
focusing behavior depends on the strength of positional disorders. The theoretical results compare 
favorably with the experimental observations, reported by Cervera, et al. (Phys. Rev. Lett. 88, 
023902 (2002)). 



PACS numbers: 43.20.- 



43.58+z, 43.90.+V 



I. INTRODUCTION 

Photonic crystals jl], ^| (PC) are made of periodically 
modulated dielectric materials, and most sonic crystals 
H (SC) are made up of materials with periodic variation 
of material compositions. Both the photonic and sonic 
crystals have been studied both intensively and exten- 
sively. It has been suggested that PCs may be useful 
for various applications such as antennas optical fil- 
ters [||, lasers ||, prisms J^, high-Q resonant cavities ||, 
wave-guides ||, mirrors JlO(] and second harmonic gener- 
ations And SCs could be used as sound shields and 
acoustic filters |J, 111 13 H 11= @> 11 H> @ • These 
applications mostly rely on the existence of photonic and 
sonic bandgaps, and a majority of them is not concerned 
with the linear dispersion region well below the first gap. 
Thus most of earlier studies were focused on the forma- 
tion of bandgaps and the propagation inhibition of waves. 

Recently, the interest in the low frequency region, 
where the dispersion relation is linear, has just started. 
Since the wavelength in this region is very large com- 
pared to the lattice constant, the wave sees the media 
as if it were homogeneous, in analogy with wave propa- 
gation in normal media. Consequently, it has been pro- 
posed j2lj that the optical lenses could be developed by 
use of PCs. However, no physical realization has been 
made so far. Along the same line of the thought, it was 
suggested that SCs may also be used to build acoustic 
lenses to converge the acoustic waves. A necessary con- 
dition to be satisfied for constructing an acoustic lens 
is that the acoustic impedance contrast between the SC 
and the air should not be large; otherwise acoustic waves 
will be mostly reflected. Once this condition is satisfied, 
the converging lens can be either convex or concave de- 
pending on whether the sound speed in the SC is smaller 
or greater than that in the air. Recently, it has been 
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experimentally found |22[ that the sound speed in two 
dimensional regular arrays of rigid cylinders in air is less 
than that in air and the impedance contrast between the 
structure and air is not large. The authors presented 
the physical realization of two refractive devices, namely, 
a Febry-Perot interferometer and an acoustic convergent 
lens. 

In this paper we carry out numerical simulation on the 
focusing of acoustic waves by SCs. We wish to theoret- 
ically analyze the experimental results reported in p2[ . 
Since the multiple scattering technique has been success- 
fully applied earlier fl4| to reproduce some experimental 
results on acoustic propagation and scattering in SCs, we 
will use this technique to study in detail about the focus- 
ing effect of acoustic waves by sonic crystals. There are 
several parameters associated with the SCs, such as the 
shapes and the arrangements of the lattice, the acous- 
tic frequencies, and the filling factors etc. We will show 
how the focusing is controlled by these parameters. In 
addition, the sensitivity of focusing effects to positional 
disorders will also be studied. The paper is organized as 
follows. The formalism is presented in Sec. II, followed 
by the results and discussions in Sec. III. Here we will 
reproduce theoretically the experimental observations in 
p2j. A summary concludes the paper in Sec. IV. 



II. FORMALISM 

Consider N straight cylinders located at fl with i = 
1,2, ...N to form either a completely random or regular 
array. An acoustic line source transmitting monochro- 
matic waves is placed at rj. The scattered wave from 
each cylinder is responsible to the total incident wave 
composed of the direct wave from the source and the 
multiply scattered waves from other cylinders. The final 
wave reaching a receiver located at ry is the sum of the 
direct wave from the source and the scattered waves from 
all the cylinders. Such a scattering could be solved ex- 
actly, following Twersky p3[. The essential procedure is 
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presented below and the detail could be found in Ref. [[15| . 

The scattered wave from the jth cylinder can be writ- 
ten as 

oo 

p s (?,f j ) = J2 ™AiH£\k\r-r j \)e in 'f>^, (1) 

n— — oo 

where k is the wave number of the medium, Hn is the 
nth-order Hankel function of first kind, and 4>p-f". is the 
azimuthal angle of the vector r—fj relative to the positive 
x axis. 

The total incident wave around the zth scatterer, 
p\ n {f), which is a superposition of the direct incident 
wave from the source and the scattered waves from all 
other scatterers, can be expressed in terms of the Bcssel 
function (of first kind) as 

oo 

PLW= E BMk\r-n\)e m ^. (2) 

n— — oo 

Using the addition theorem for Bessel function, the scat- 
tered waves p s (f,fj) for each j =/= i can be expressed in 
terms of the modes with respect to the ith scatterer and 
thus one leads to 

N 

^=4+ E (3) 

with 

Sf = iirH^ (k\ - il^fil) (4) 

and 

oo 

CiC= E i7rAfffW B (fe|r1-r7|)exp[ia-n)^i_ rl ]. (5) 

/ — — oo 

Matching the usual boundary conditions for each cylin- 
ders, one obtains: 

B* n = i7rP„4 (6) 

where T l n are transfer matrix relating the acoustic prop- 
erties of the scatterers and the media (see Ref. for the 
expression). All coefficients A l n may easily be obtained 
by solving the Eq. (^) along with Eq. (|6|) and hence the 
total wave at any desired point outside the cylinders may 
be obtained as 

N oo 

p(r)=Po(r)+E E ^KH^(k\r-r t \)e m ^. (7) 

i— 1 n— — oo 

The acoustic intensity (field) is the square module of the 
transmitted wave where the normalized transmission is 
given as T = p/po- It is worth to mention that the to- 
tal wave expressed by Eq. (]?]) incorporates all orders of 
multiple scattering and the formalism is valid for both 
regular and random arrangements of the scatterers. 



III. RESULTS AND DISCUSSIONS 

In the following computation, we consider regular ar- 
rays of rigid cylinders in air. The system we use is from 
and has been detailed in j2^] . There are N uniform rigid 
cylinders of radius a. The area occupied by the cylinders 
per unit area is (3. Therefore, the lattice constant can 

be expressed as d — y (2ira 2 /y/3f3) for a triangular lat- 
tice arrangement of the cylinders. For the corresponding 
random arrangement of the cylinders, d represents the 
average distance between the nearest neighbor cylinders. 
Unless otherwise mentioned, the lengths will be scaled by 
d and the frequency will be scaled as ka, to make them 
dimensionless. We will mainly consider the triangular 
lattice and use the parameters from |^2| for the filling 
factor, and the radii of the cylinders. 




Band Structure 



FIG. 1: Band structure for the two dimensional triangular 
lattice of the rigid cylinders in air. 

Figure [j] reproduces the band structure for the two 
dimensional triangular lattice from (^2|. The filling frac- 
tion is 0.2024, the radius of each cylinder is 1.5cm, and 
thus the lattice constant is « 6.35cm. We see that the 
dispersion in the low frequency region is linear for both 
the directions. Our calculations will be confined in the 
linear low frequency region only. 

In accordance with the experimental setup described 
in the conceptual setup of the source and the lat- 
tice arrangement of the cylinders for the computation 
are shown in Fig. |^(a). The line source is denoted by 
S. Here the triangular arrangement in the experiment is 
taken, and the lenticular shape is identical to that used in 
the experiment [Q. The lattice constant for the lattice 
is ~ 6.35cm, the same as considered in the experiment. 
The source is placed at lOOd away from the left side of 
the lattice; this distance is far enough to ensure that the 
incident wave on the lattice behaves as a plane wave. The 
wave with frequency of ka — 0.486 (~ 1700iJz) is inci- 
dent on the triangular lattice of lenticular shape from the 
left, and the transmitted intensity (|T| 2 ) is calculated on 
the right side of the lattice. The two dimensional spa- 
tial distribution of the transmitted intensity is shown in 
Fig. @(b). Here the x axis is a horizontal line towards 
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right, and the y axis is placed vertically upward. The 
origin is at a distance twice the radius of the cylinders 
to the right most cylinder, so to avoid the diverge when 
plotting the spatial distribution. The same arrangement 
of the coordinates is also used in the following figures, 
except for Fig. We see that the focusing of the trans- 
mitted wave is evident, and is in agreement with the 
experimental observation ]2^| . Since the maximum in- 
tensity point, which we assume as the focal point, is not 
clear from Fig. ||(b), we have computed the variation of 
the field intensity along the x axis with y = 0, and also 
along the y axis with x = 20. The results are shown in 
Figs. |^(c) and ||(d). Fig. ||(c) shows that there is a peak 
at x ~ 14. The intensity variation along the y axis shows 
a significant peak at y = 0. Here we observe that the in- 
tensity is better confined along the y axis than along the 
x axis. The theoretical results are in a good agreement 
with the experimental observation. 




Y axis 



FIG. 2: (a) The line source denoted by S and the arrangement 
of cylinders; (b) two dimensional spatial distribution of the 
transmitted intensity (|T| 2 ) on the right side of the crystal; 
(c) the variation of the intensity in (b) along the x axis at 
j/=0 and (d) the variation of the intensity in (b) along the y 
axis at x=20. 



In view that our numerical calculation can reproduce 
the experimentally observed focusing effect of sound, we 
are encouraged to further explore possible focusing ef- 
fects with different shapes of SCs. Here we consider the 
lenticular shape shown in Fig. ||(a). The size of the lens 
along the y direction is same as that for Fig. ||, whereas 
the size along the x axis is larger. Thus, the edges of the 
lattice are much smoother than that in the above lens. 
Here the source is again kept at lOOd away from the left 
side of the lens, and other parameters are the same as 
used above. The transmitted intensity calculated on the 
right side of the lens is shown in Fig. ||(b). We see that 
the focusing is apparently better than that in the above 
SC of lenticular shape. We also plot the intensity vari- 
ation along both the x axis with y = and the y axis 
with x — 2. The results are shown in Figs. ||(c) and 
||(d) respectively. Fig. [|(c) clearly shows that the inten- 
sity is peaked at about x = 2 and y = 0, then decreases 



more rapidly compared to that in Fig. 0(c) . Fig. 0(d) also 
shows a narrower peak when compared with Fig.^d). In 
this case, the focal length is shorter and the intensity at 
the focal point is much stronger than those in the above 
case. Again, the intensity region is more extended along 
the x axis compared to that along the y axis. 




X axis 



Y axis 



FIG. 3: (a) The arrangement of the cylinders; (b) the two 
dimensional spatial distribution of the transmitted intensity 
(|T| 2 ) on the right side of the lens; (c) the variation of the 
intensity in (b) along the x axis at y=0 and (d) the variation 
of the intensity in (b) along the y axis at x=2. 

For comparison, we have also studied the transmission 
through a rectangular lens shown in Fig. ||(a). The re- 
sults are depicted in Fig. [|(b). In this case, no unique 
focusing is observed. The transmitted intensity is more 
or less spreading all over the xy plane. This can be better 
seen from Figs. |](c) and|^(d), where the intensity varia- 
tions are plotted along the x axis and y axis respectively. 
These results also agree with the experimental observa- 
tion 01. 




FIG. 4: (a) The arrangement of the cylinders; (b) the two 
dimensional spatial distribution of the transmitted intensity 
(|T| 2 ) on the right side of the crystal; (c) the variation of the 
intensity in (b) along the x axis at y = 0, and (d) the variation 
of the intensity in (b) along the y axis at x = 5. 

We also performed the similar calculations for the lat- 
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tices with the regular square arrangement of the cylinders 
and found that the results are similar. 

Figure ||(a) shows the variation of the focal length, 
i. c. the distance between the focal point and the center 
of the lattice, as a function of the filling factor. We see 
that the focal length (£ in cm) decreases rapidly with the 
increasing filling factor (j3) up to f3 = 0.2, then tends 
to saturate for higher filling factors. This decrease of the 
focal length with the filling factor is reasonable as the the 
effective sound speed in the SC decreases with increase 
of the filling factor, to be discussed later. The variation 
of the intensity at the focal point (I max ) as a function 
of filling factor is also shown in Fig. ||(b). Therefore, 
in addition to the shapes, the filling factor is another 
parameter that effectively controls the focal point and 
focusing intensity. 
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FIG. 5: (a) The variation of the focal length as a function of 
filling fraction and (b) the variation of Maximum inten- 
sity, Imax (the intensity at the focal point) as a function of 
the filling fraction. 

From the above study on lattices of different shapes, 
we conclude that lattice shape is very crucial to obtain 
the focusing effect. The quality of focusing and the focal 
length could be adjusted by varying the shape of the 
crystal or the filling fraction. 

In practice, it is useful to know the robustness of the 
focusing against possible disorders in the system. Here 
we consider the positional disorder of the cylinders. Since 
the wavelength is larger than the lattice constant in the 
present cases, one might intuitively conclude that the po- 
sitional disorder has no effects. This is untrue. We have 
performed the calculation for three cases: the weak, mod- 
erate, and strong positional disorders. For the weak and 
moderate disorders, every cylinder can randomly take a 
position within a circle with radius of 1.4 and 15 percent 
of the lattice constant about its regular position, while 
for the strong disorder, the cylinders may take positions 
randomly within the boundary of the regular crystal as 
long as no two cylinders overlap with each other, i. e. 
completely random. Figs. |^(a) and (b) show the two di- 
mensional spatial distribution of the transmitted inten- 
sities for the weak and moderate disorders. It is clear 
that the focusing effect is hardly affected by the presence 
of weak disorder. Fig. ^(c) shows the distribution of the 
transmitted field for the strong disorder case. Here the 
the intensity is highly suppressed all over the xy plane 
and also the intensity at various places in the xy plane 
becomes comparable to that at the maximum intensity 
region. Therefore, we conclude that slight or mild de- 



viations from perfectly regular arrangements when de- 
signing acoustic lenses will not affect the focusing dras- 
tically; therefore it will be easier to build the acoustic 
lens. However, the strong disorder does affect the focus- 
ing significantly. It is interesting to note that although 
the focusing nature is greatly reduced in the presence of 
strong disorders, the spatial pattern of the scattered in- 
tensity remains almost unchanged. This is because at the 
long wavelength region, the waves will not be able to dis- 
cern the detailed structure of the crystal, and therefore 
the scattering pattern is mainly controlled by the shape 
of the crystal. The reduction in the transmitted intensity 
indicates that the waves are reflected backward. This is 
confirmed below. 




X axis 

FIG. 6: Two dimensional spatial distribution of the transmit- 
ted intensity (|T| 2 ) on the right side of the lattice with (a) 
weak, (b) moderate, and (c) strong positional disorders. The 
shape of the lattice and the filling factor used here are the 
same as that depicted for Fig. ^. 

To understand such drastic difference for the focusing 
between the regular lattice and the disordered lattice we 
plot the spatial distribution of back scattered intensities 
for the regular and the disordered lattices and they are 
shown in Fig. 0(a) and 0(b) respectively. Here the x axis 
is a horizontal line towards right, and the y axis is placed 
vertically upward. The origin is at a distance twice the 
radius of the cylinders to the most left cylinder. We ob- 
serve that the distribution is symmetric for the regular 
case, but asymmetric for the disordered case. Compared 
to the forward scattering pattern for the disordered case 
in Fig. 0, the backscattering is more sensitive to the ran- 
dom configuration. This is understandable since when 
waves are backscattered, the travelling path will be twice 
that when going forward, thereby causing more interfer- 
ence between scattered waves. The intensity for the dis- 
ordered case appears to be stronger than in the regular 
case. The quantitative difference is further illustrated by 
Figs. 0(c) and 0(d) where the variations of back scat- 
tered intensities for both the regular and disorder cases 
are plotted along the x axis with y — and the y axis 
with x=0 respectively. While the solid curve represents 
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the variation of back scattered intensities for the regular 
case, the dashed curve for the disordered case. Here, it 
is clear that the intensity strength for the disorder case 
is much higher. In fact, the ratio of the maximum back 
scattered intensity between the two cases can be as high 
as 15. In other words, in the disordered configuration, 
the transmission is highly prohibited, in agreement with 
the earlier observation [jfi]. In this case, the lattice does 
not appear to be transparent to acoustic waves. 



effective medium is solved as 



Pe 



1-/3 



(9) 



For a given pressure, the sound speed in an air filled 
medium is proportional to the inverse of the square root 

of the density, i. e. c ~ ■\J~^, thus we obtain the sound 

speed of the sonic crystal 
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FIG. 7: Two dimensional spatial distributions of the back 
scattered intensities from (a) the regular lattice, and (b) one 
completely random configuration of the lattice; no configura- 
tion averaging is taken. The shape of the crystals for both 
cases is the same as that in Fig. 0. 



(10) 



Obviously, by 1 — p « 1/(1 + /?) for small P, the result 
reduces to that in (22) for small filling factors. In Fig. ||, 
we show the comparison among the experimental data, 
and the theoretical results from the present model and 
that in^2|. Here we clearly see that present modified 
model reproduces the experimental data remarkably well, 
better than the original model. 




Finally we note on the acoustic refraction from the 
sonic crystal formed by the rigid cylinders in air. The 
focusing effect shown above is essentially, as already ar- 
gued by the authors in [53 , due to the reduction of the 
sound velocity within the SC. A simple model has been 
proposed by the authors to explain this reduction. Upon 
examination, we found that the results from the model 
agrees with the experimental data rather well in the low 
filling regime. A deviation between the theoretical and 
experimental data is obviously seen, referring to Fig. 2 in 
p2[ , for high filling factors. Here we would like to pro- 
pose a modified model to account for the experimental 
data. At the frequency considered in the experiment Q 
and the simulation above, the wavelength is much bigger 
than the lattice constant. Therefore the sonic crystal may 
be regarded as an effective air-filled medium, i. e. the air 
around each cylinder is regarded as an effective air which 
takes into account of the presence of other cylinders. The 
presence of the rigid cylinders is actually to reduce the 
occupation of the air per unite area without the cylin- 
ders, and is replaced by the effective air medium. This 
scheme may be called the self-consistent procedure. The 
effective mass density can be then estimated from 



Pe = Pair + PPe, 



(8) 



where p a i r and p e are the densities of the air and the 
effective medium respectively. Thus the density of the 



FIG. 8: The variation of sound speed in the SC with the 
filling fraction: The circles with error bars represent the ex- 
perimental values, the dashed curve represents the theoretical 
estimate in [123] , and the solid curve represents our estimate. 



IV. SUMMARY 

Through rigorous calculations we have successfully re- 
produced the focusing of acoustic waves by the sonic crys- 
tal made of rigid cylinders in air, observed by the recent 
experiment. We have examined the focusing effect for 
crystals of different shapes. It has been shown that the 
shape of the crystal plays a crucial role on the quality 
of focusing. The dependance of the focal length on the 
filling factor is also examined. It is found that the focus- 
ing becomes brighter and the focal length decreases as 
the filling factor increases. The robustness of the focus- 
ing against the presence of disorders is also investigated. 
The results show that the robustness of the focusing be- 
havior depends on the strength of disorders. Finally, a 
model describing the sound speed in the SC made from 
rigid cylinders arranged periodically in air as a function 
of filling fraction is proposed. An excellent agreement 
is found when the results from the model are compared 
with the experimental data. 
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